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Abstract 

We obtain the equations that define the equilibrium of a homogeneous relativistic 
gas of neutrons, protons and electrons in a constant magnetic field as applied to the 
conditions that probably occur near the center of neutron stars. We compute the 
relative densities of the particles at equilibrium and the Fermi momentum of electrons 
in the strong magnetic field as function of the density of neutrons and the magnetic 
field induction. Novel features are revealed as to the ratio of the number of protons to 
the number of neutrons at equilibrium in the presence of large magnetic fields. 
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Significant interest in the behaviors of relativistic electrons in a constant strong magnetic 
field was greatly stimulated by the discovery of giant magnetic fields at the surface of neutron 
stars 0, 0] . Observations of hard X-rays from pulsar Hercules X-l allowed one to estimate the 
magnetic field induction at the pulsar surface to be of the order of 10 13 G. Such a magnetic 
field "frozen" in a neutron star must become stronger and stronger reaching the ultrastrong 
value of the order of 10 18 G in the central part of the neutron star ||. Furthermore, the 
fields of the order of 10 23 G [|J and, even possibly 10 33 G || may exist at the electroweak 
phase transition. It is worthy of note that the magnetic field in a neutron star may be 
considered as a macroscopically uniform field with respect to the characteristic scales, such 
as the Compton wavelengths, of constituent particles of the neutron star. However, the field 
is extremely nonuniform at the electroweak phase transition. In any case, such a strong field 
will affect in an essential way the behaviors of charged particles in the star. 

In this paper, we discuss the effects of strong constant magnetic fields on the chemical 
equilibrium of a degenerate gas of neutrons, protons, and electrons. We assume in what 
follows that the gas is spatially homogeneous due to the homogeneity of the magnetic field. 

In a very strong magnetic field with induction B higher than Bq = m 2 c 3 /eh = 4.41 • 10 13 
G (m and e are the electron mass and charge, respectively), electrons in the lowest Landau 
energy level "move" in the plane transverse to the magnetic field direction in the region with 
characteristic dimensions of the order of %{Bq/ B) 1 ^ 2 jmc (see, e.g. ||). We note that each 
electron energy level is degenerate with an infinite multiplicity . However, for the case under 
consideration the number of quantum electron states is restricted by both the value of the 
magnetic field induction and the size of the region in which the magnetic field acts. 

The first physical quantity which we must define is the chemical potential /i e , or the 
Fermi energy (if the gas temperature is equal to zero), of relativistic electrons as a function 
of magnetic field induction (see also [[(J). The electron number density n e in the presence of 
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a magnetic field B = (0, 0, B) is defined by 

eB x , r°° dp 
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where _E n = (p 2 c 2 + m 2 c 4 + 2eBhcn) 1 ^ 2 is the energy spectrum of a relativistic electron M; 
n = 0, 1,2, . . . enumerates the Landau levels, ( is the quantum number characterizing the 
spin projection of an electron, p is the electron momentum component parallel to B, and 
9 is the temperature. In what follows we shall express values of /i e in dimensionless units: 
(jl = He/mc 2 . Since the chemical potential /i of an electron gas at 9 — and B = is related 
to n e by 

(/i 2 - 1) 3/2 = n e /n , n = m 3 c 3 /3n 2 h 3 , (2) 

eq. ([I]) may be considered as a relation between /i, /x and # at given B. By integrating with 
respect to p in (0), taking allowance for -E max = /i e at # = 0, it is easy to obtain: 
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The value n max is determined by the condition n max < (fi 2 — 1)B /2B. 
It is seen from that, if the following condition is satisfied: 

(V2n c /3n ) (B /B) 3/2 < 1 , (4) 

then only the lowest (ground) level n = remains populated. For a given n e , we can find 
from (|^) the relation 

( /i 2 _l)V 2 = (2B Q /3B)(n /no). (5) 

Let us now consider the conditions for chemical equilibrium of a degenerate gas of protons 
(p), neutrons (n) and electrons (e) in the presence of large magnetic fields. First we shall 
consider fields B in the range 

B* > B > B , (6) 



where Bq = m 2 c 3 /eh = 3.4- 10 6 Bo, and m p is the proton mass. Magnetic field B in this range 
is "strong" for the electron but "weak" for the proton. Hence motion of protons in such a 
magnetic field can be considered as quasiclassical because the spacing between Landau levels 
for the proton under these conditions is very small. We assume further that the condition 
(f§) is also satisfied. Then only the lowest energy level will be occupied by electrons. We 
also suppose that the temperature 9 is equal to zero, since for a typical 100-year old neutron 
star, its temperature is estimated to be 10 8 K (about 10 keV), which can be considered cold 
as compared to the Fermi energy (about 1000 MeV) of the degenerate relativistic neutrons 

We are interested in the reactions in which the total density of baryons = n p + n n is 
conserved and the electroneutrality condition of a gas n p = n c is satisfied . These processes 



are called the direct URCA processes |9], [10], [LI], [T!|. Since n b is conserved, the total energy 
density e depends only on n n : 

e = n- 2 h- 3 f* p 2 (p 2 c 2 + mYfl 2 d V + tT 2 ^ 3 f * p 2 (p 2 c 2 + m 2 J) x l 2 dp + e e (p c F , B), (7) 
Jo Jo 

where 

Pf = (n n /n 0n ) 1/3 m n c, = (n p /n 0n ) 1/3 m D c, 
p e F = {2n p /3n ){B /B)mc, n 0n = m 3 n c 3 /37r 2 h 3 , (8) 

and e e is the energy density of the electrons. 

At chemical equilibrium de/dn n = 0, or /x n = /x p + /i e , where /i n ,/i p and /x e are the 
chemical potentials of neutrons, protons and electrons respectively. The chemical potential 
of a degenerate electron gas in the presence of the magnetic field B under the conditions 
assumed is /i e given by formula @. The quantity pp plays the role of the Fermi momentum. 
More precisely, p F here is the "Fermi momentum projection" in the magnetic field direction 
for electrons. 



The chemical potentials of relativistic nucleons are given by 
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for B satisfying (|J). The proton number density n p as a function of n n at chemical equilibrium 
can be found from the following equation: 
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Writing inequality (f|) in the form 



n p /n 0n < (3/v^) (m/m n ) 3 (B/B f 2 « 3.4 ■ 1(T 10 (B/B ) 
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we see that, for £?q 3> -B 3> -Bo, the inequality n p /n Qn 1 must be satisfied. 

We now turn to discuss equilibrium condition for the case of ultrastrong magnetic fields 
(with B > Bq) and very high densities (with n p , n n > no n ). The existence of such ultrastrong 
magnetic fields near the center of neutron stars is allowed, or, at least, is not forbidden in 
principle. In this case formula (|T0"D needs to be modified because now the protons, just as the 
electrons, must occupy their lowest Landau energy levels. The limit of the proton density 
under which all protons, and consequently all electrons, populate only their respective ground 
levels is again given by (JTTJ) . But now one could have n p /n Qn > 1. As in the case of electrons, 
the Fermi energy for protons is given by 
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Eq. fllOD is then replaced by 
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Numerical solutions of eq. (|10D and (|T3|) subject to the constraint (]TT|) are given in Fig. 1 
(B < 10 5 B Q for eq. ([T0|) , and B > 10 7 -Bo for (|T3|) ). These curves represent the normalized 
proton density number n p /rio n as a function of the normalized neutron density number 
u n J n 0n at equilibrium at various values of magnetic field induction. The dotted curve gives 
the corresponding values in the absence of magnetic fields. The normalized Fermi momenta 
of electrons pp / mc as a function of the normalized neutron density number for these values 
of magnetic field are given in Fig. 2. 

One can see from Fig. 1 that results obtained here differ from analogous results with zero 
magnetic field. Such a difference takes place at various values of magnetic field induction 
for some ranges both of low and high particles densities. Under the conditions assumed in 
this paper, for a given value of n n /n 0n , the values of n p /n 0n in the presence of finite B's are 
for the most part higher than the corresponding value when the field is absent, until the 
proton density reaches a value (n p /n 0n ) cross = (3B /2B ) 3 ^ 2 (m/m n ) 3 , which is the point of 
intersection between the curves with B ^ and B = 0. The value of (n p /rion)cross is slightly 
less than the upper limit given in ( pTTj) for a given B. Also, in the presence of ultrastrong 
magnetic field, there appear values of densities for which n p > n n . These ranges of particle 
densities are of great interest since in zero magnetic field the ratio n p /n n is always less than 
unity, with a maximum equals 1/8 [10|. The functional form of (|1^) implies that n p /n 0n is 
directly proportional to B. This is evident in Fig. 1 from the fact that curves corresponding 
to B > 10 7 B are parallel to one another, and are seperated by equal distance. 

One may understand the behavior of these curves as follows. Eq. (H) shows that higher 
magnetic field B tends to lower the Fermi momentum p| of the electrons, and hence the 
chemical potential of electrons (chemical potential of protons as well, when B > Bq, see 
(|12"D). To maintain chemical equilibrium among the particles at fixed value of n n /no n , the 
chemical potentials of electrons and protons have to be raised so that eq. fllOl) or (|T3|) is still 
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satisfied. This is achieved through the increase in the density of electrons, and hence the 
density of protons by neutrality condition assumed here. The shape of the curves, which 
consists mainly of segments of straight lines in the log- log plot, indicates power- law behavior 
between the quantities n p /n 0n and n n /n Qn . This can also be easily understood. Consider first 
the situation in which B > Bq. In the high density region where n n /no n ^> 1, the densities 
of electrons and protons are high enough that their Fermi momenta satisfy p F 3> mc and 
Pp m p c. Hence only the second term within each square-root in flUf ) dominates. After 
some algebra we get 



n 
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This shows that the straight lines on the upper right part of Fig.l have slope 1/3. For 
the region where n n /riQ n <C 1, the densities of electrons and protons are lower such that 
p F -C m p c, but still pp 3> mc. Eq.(|I3|) then reduces approximately to 



m n 

2 \n 0n 

From this one has, after setting m n w m p , 
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The slopes of the corresponding line segments in Fig.l are 2/3. We note here that (|i~4D 
and (|I~6"D have the same coefficients in front of the factor n n /n Qn , which mean the two line 
segments in Fig.l have the same intercept. Similar consideration, when applied to ( [10| ) for 



the case Bq 3> B 3> B , leads also to the power law fllBl). When B = 0, we get the power 
laws — = \— and — = |(— ) 2 for high and low neutron density, respectively. 

Fig. 2(a) shows that, for a given neutron density, the Fermi momentum of electrons in the 
presence of field B < Bq decreases as B increases. It is not higher than the corresponding 
value when B = until (n p /n 0n ) cross is reached. On the other hand, one sees from Fig. 2(b) 



that electron Fermi momenta in the presence of ultrastrong magnetic fields are independent of 
the magnetic field strength. Furthermore, they equal the corresponding value in the absence 
of magnetic field. The former situation observed in Fig. 2(b) is just the manifestation of 
the fact that, for a fixed n n /no a value, n p /rio n is directly proportional to B in this case, 
as mentioned before, and the fact that p F depends only on their ratio. To understand the 
latter result, we simply observe that the equilibrium equation fllcf) , expressed in terms of the 
Fermi momenta of the nucleons, and the electrons, is identical to the corresponding equation 
without magnetic field. This is because p F has the same form in both cases, and while the 
corresponding expressions of p F and p F differ in the two cases, the relation p F = p F holds 
well. So for a fixed neutron density (i.e. fixed p F ), both equations give the same solution for 
the value of pp, solved by different n p /no n at different B when B ^ 0, of course. 

The results obtained here may serve to give us hints as to the chemical composition of the 
central region of neutron stars, if we suppose that this region is composed of "cold" degenerate 
neutrons, protons and electrons in the superstrong magnetic field. The assumption was also 
made that all the electrons occupied only the lowest energy level when B ^> Bq. If the gas 
temperature in this region is not equal to zero, then we must introduce definite restrictions 
not only on the quantity n e but also on the gas temperature. There are several aspects need 
to be considered. Firstly, it is clear that at 9 ^ the electrons can occupy the higher energy 
levels of transverse (to the magnetic induction vector) motion of electrons due to thermal 
excitations. Secondly, the electron-positron pair production due to electron collisions in 
the ultradense hot medium can become noticeable. Finally, the chemical potential of the 
electrons depends on the gas temperature. 

It is easy to see that the electrons will occupy the lowest energy level of transverse motion 
at B > B if 




(17) 



8 



where ks is the Boltzmann constant. It follows from this inequality that, even at a very high 
temperature of the order m , electrons are still in the lowest energy state. One must also 
estimate the effect of electron-positron pair production by electron collisions in the presence 
of a magnetic field. As a result of electron-positron pair production, the equality n p = n e will 
be violated but the electroneutrality condition as a whole must be conserved. We estimate 
that, even when B 3> Bq, the densities of electrons and positrons produced will be much 
less than the initial density of electrons n e if the temperature satisfies the strong inequality 
kBO/mc 2 <C (n e Bo/3no n B In 2) (m n /m) 3 . Thus the thermal effects of the electron-positron 
pair production (even when ksO is of order mc 2 ) will not substantially affect the chemical 
equilibrium of neutrons, protons and electrons over a sufficiently wide interval of values n e . 
We have also made an estimate of the chemical potential of electrons which occupy the 
lowest Landau energy level as a function of temperature and found that, if inequality fllTD 
is satisfied, then the chemical potential in the form (|5|) can still be employed. 

Finally we remark that in this paper we consider only the case of ideal gas model, as 



are usually done in the general analyses of basic properties of compact stars || [10| . We 
see, however, that even such simple model gives nontrivial and unexpected results. So a 
careful analysis of more realistic models in the presence of strong magnetic fields, and at 
finite temperatures, is extremely needed to answer the question about the actual behavior of 
such models. Many factors need be considered. For instance, the density of matter near the 



center of a neutron star is probably an inhomogeneous function of distance |TT], [L2[], and the 
problem of equilibrium in this case is significantly more complicated. Also, one must consider 
the effect on the chemical equilibrium by the coulomb interaction and the quantum exchange 
effect among the particles. These consideration will be reserved for future investigations. 
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Figures Captions 



Fig. 1. Plot of log(n p /n 0n ) verses log(n n /n 0n ) (continuous curves) for values of B/B Q 
at (from bottom to top): 10 2 , 10 3 , 10 4 , 10 5 , 10 7 , 10 9 , 10 11 and 10 13 . The dotted line 
corresponds to the curve with B = 0. Equation log(n p /n 0n ) = log(n n /n 0n ) is indicated by 
the dashed line. 

Fig. 2 (a). Plot of the normalized Fermi momenta of electrons pp/mc verses n n /n Qn 
for B/Bq = 10 4 (upper curve) and 10 5 (lower curve). Dotted curve corresponds to B — 0. 

Fig. 2 (b). Plot of the normalized Fermi momenta of electrons for p F /mc verses n n /n 0n 
for B/B = 0, 10 7 , 10 9 , 10 11 and 10 13 . 
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